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In situ Evidence of Breaking the Ion frozen-in Condition via the 
Non-gyrotropic Pressure Effect in Magnetic Reconnection 
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For magnetic reconnection to proceed, the frozen-in con¬ 
dition for both ion fluid and electron flnid in a localized 
diffusion region must be violated by inertial effects, thermal 
pressure effects, or inter-species collisions. It has been un¬ 
clear which underlying effects unfreeze ion fluid in the diffu¬ 
sion region. By analyzing in-situ THEMIS spacecraft mea¬ 
surements at the dayside magnetopause, we present clear ev¬ 
idence that the off-diagonal components of the ion pressure 
tensor is mainly responsible for breaking the ion frozen-in 
condition in reconnection. The off-diagonal pressure tensor, 
which corresponds to a nongyrotropic pressure effect, is a 
fluid manifestation of ion demagnetization in the diffusion 
region. From the perspective of the ion momentum equa¬ 
tion, the reported non-gyrotropic ion pressure tensor is a 
fundamental aspect in specifying the reconnection electric 
field that controls how quickly reconnection proceeds. 


1. Introduction 

Magnetic reconnection is considered to drive global-scale 
dynamics in Earth’s magnetosphere [Angelopoulos et al, 
2008, 2013], solar flares [Parker, 1957; Petschek, 1964], lab¬ 
oratory plasmas and astrophysical systems. Magnetic re¬ 
connection changes the magnetic field topology and releases 
magnetic energy into particle energy in plasmas. As recon¬ 
nection occurs, magnetic field lines ‘break’ and ‘reconnect’ 
at the X-line. Outside the diffusion region, plasma motions 
are frozen to magnetic held lines that behaves like elastic 
strings. Within the diffusion region, the violation of plasma 
frozen-in allows magnetic field lines to disconnect and re¬ 
connect. 

The name ‘diffusion region’ originates from the early 
Sweet-Parker reconnection model in which oppositely di¬ 
rected magnetic field lines (B) in the current sheet diffuse 
into the plasma (cltB ~ V^B) as plasmas are demagnetized 
by inter-species collisions [Parker, 1957; Petschek, 1964]. 
The term ‘diffusion region’ has been extended to include any 
physical process that can violate the frozen-in condition [ Vo- 
syliunas, 1975; Hesse et al., 1999]. The physics in the limited 
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diffusion region is of high importance in magnetic reconnec¬ 
tion. The global-scale evolution of the magnetic topology 
relies on the reconfiguration of magnetic field lines in the 
diffusion region. The mechanisms responsible for breaking 
the frozen-in condition specify the reconnection electric field 
that controls how quickly reconnection proceeds. 

For reconnection to proceed, the frozen-in condition for 
both ion and electron must be violated in the diffusion re¬ 
gion. In case that only one species is unfreezed from the field 
line, one can trace the frozen magnetic field lines tied to the 
other magnetized species. The exact condition for unfreez¬ 
ing ion fluid is a non-zero curl of E + Vi x B based on the 
frozen-in theorem. According to the ion fluid momentum 
equation without approximation, 

E-l-ViXB = {l/niqi)\/-Pi+{mi/qi)d\ri/dt+{mi/qi)uie(vi—Ve), 

( 1 ) 

where E is the electric field. Pi is the ion pressure tensor, 

Vi (ve) is the ion (electron) flow velocity, and Vie is the ef¬ 
fective ion-electron collision frequency. As evident from the 
ion momentum equation (1), the ion frozen-in condition can 
be violated by three non-ideal effects, i.e., an anisotropic 
ion pressure tensor, an ion inertial (acceleration) effect or 
inter-species collisions manifested as friction. These non¬ 
deal effects are the fluid manifestation of the kinetic effects 
leading to ion demagnetization in reconnection. The expres¬ 
sion (1) and arguments for ions applies equally to electrons 
after a suitable change of charge, mass and sub-indexes. 

The main non-ideal effects for ion and electron fluid in 
reconnection is still an open question. Reconnection mod¬ 
els are categorized according to which non-ideal process 
is thought to be dominant. In the resistive reconnection 
model, inter-species collision is the dominant effect. The col¬ 
lisions can be either binary or anomalous (induced by wave- 
particle interactions that are not clearly understood). In 
collisionless reconnection models where effective ion-electron 
collisions are infrequent, the non-ideal effects have to be 
anisotropic pressure effects and/or inertial effects. As plas¬ 
mas are convected toward the reconnection site, ions first 
become demagnetized at the characteristic scale length of 
ion motions while electrons are still magnetized [Sonnerup, 
1979]. Electrons are expected to be demagnetized at a 
smaller spatial scale. 

The reconnection electric held, defined in the co-moving 
frame of the current layer, controls how quickly (by E- J) 
magnetic energy is converted to particle energy in recon¬ 
nection. The reconnection electric held is the same in the 
ion momentum equation (Eq.l) and the electron momen¬ 
tum equation. From the perspective of the ions (electrons), 
the reconnection electric field is the sum of the non-ideal 
electric held supported by ion (electron) non-ideal effects 
and the simple ion (electron) convection electric field. Al¬ 
though demagnetization is more difhcult for electrons, the 
reconnection electric held is equally explainable in terms of 
ions or electrons near the X-line [Hesse et al., 2011; Cai and 
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Lee, 1997]. In this sense, ion demagnetization and electron 
demagnetization are two different but equally fundamental 
aspects in specifying the reconnection electric field. 

In collisionless reconnection, ion decoupling from the 
magnetic field lines is considered to produce a Hall effect 
\Sonnerup, 1979; Nagai et ai, 2001; 0ieroset et al., 2001]. 
But what non-ideal effect is responsible for unfreezing ion 
fluid in the first place has not been clear. Intuitively, it 
is tempting to attribute ion decoupling to the Hall current 
term in the generalized Ohm’s law, E-fVi x B = (l/ng)J X B. 
But this relation is a trivial equivalence of the electron 
frozen-in condition (E-|-VeXB=0). Ion information is can¬ 
celed out from this relation and nothing specific about non¬ 
ideal effects of ions can be inferred. Studies of ion-scale 
physics in reconnection have been usually guided by the 
generalized Ohm’s law, which is essentially a combination 
of electron and ion momentum equations with some approx¬ 
imations. Although appropriate for electrons to a certain 
extent, these approximations neglect and hide important as¬ 
pects of ion dynamics. The correct approach to investigate 
ion scale physics is the full momentum equation without ap¬ 
proximation [Hesse et al, 1999, 2011; Cai and Lee, 1997]. 

Experimental clarification of the ion non-ideal effects re¬ 
quires a comparison between the non-ideal electric field 
(E-b Vi X B) and ion pressure/inertial terms. Although at¬ 
tempts have been made to compare the non-ideal reconnec¬ 
tion electric field with the divergence of the electron pressure 
tensor , the electron diffusion region was too small for the 
spacecraft to encounter [Henderson et al, 2006]. On the ion 
skin depth scale, past spacecraft observations have reported 
ion kinetic features such as non-gyrotropic ions and coun¬ 
terstreaming ions during reconnection [Hoshino et ai, 1998; 
Nagai et al, 1998; Wygant et al., 2005; Zhou et ai, 2009; 
Aunai et al., 2011]. But how these effects are linked to ion 
demagnetization and the violation of the ion frozen-in condi¬ 
tion has been elusive. Here we report observations from an 
encounter of the THEMIS spacecraft [Angelopoulos, 2008; 
Sibeck and Angelopoulos, 2008] with the diffusion region near 
the reconnection X-line at the Earth’s magnetopause. The 
comprehensive fields instrumentation and measurements of 
ion velocity distribution [Bonnell et al., 2008; McFadden 
et al., 2008] on THEMIS spacecraft provide an ideal op¬ 
portunity to address the question of ion demagnetization 
in reconnection. By comparing the non-ideal reconnection 
electric field with ion pressure/inertial terms, we for the first 
time identify the non-ideal effects corresponding to ion de¬ 
magnetization in the diffusion region. 


2. Data analysis 

On Feb 13, 2013, THEMIS spacecraft E moved into the 
magnetopause and detected a reconnection diffusion region. 
THEMIS E was at (6.1,-6.9,-0.6) in units of Re (Earth 
radii) in the Geocentric Solar Magnetospheric coordinate 
system (GSM). We adopted the boundary normal coordi¬ 
nate NML to study the magnetopause current sheet. N (- 
0.78,0.59,0.21)GSM is the boundary normal direction (out¬ 
ward) as determined by the minimal variance direction of 
magnetic fields, L=(-0.07, -0.42, 0.90) GSM is the direction 
of maximum variance of the magnetic field, and M com¬ 
pletes the right-handed orthogonal coordinate, directed out 
of plane. Figure 2 shows measurements of fields and parti¬ 
cles during the magnetopause crossing. The magnetopause 
current sheet, indicated by a change of B^ from -40 nT to 60 
nT, was observed from 23:24:40UT to 23:25:00UT. North¬ 
ward plasma flow with velocities as large as 130 km/s was 
detected from 23:24:42UT to 23:25:05UT, immediately after 


the spacecraft crossed the separatrix SI at 23:24:40UT. The 
separatrix SI is identified by a sudden increase of 1-10 keV 
electrons that originated in the magnetosphere. The elec¬ 
tron characteristics are good indication of separatrix field 
lines because their small gyroradius and large mobility along 
the magnetic field. Around 23:25:10UT, THEMIS E crossed 
the separatrix S2 and moved into the magnetosphere. The 
separatrix S2 is identified by a boundary between a broad 
spectrum of mixed electrons and a dominant magnetosphere 
electrons. Between 23:25:lOUT and 23:25:30UT, spacecraft 
was away from the current sheet and the diffusion region, as 
evidenced by the large B^ component. After 23:25:30UT, 
THEMIS E moved back into the magnetopause and observed 
the southward plasma flow. Bjv shows a variation of a few 
nT on a 20s timescale in addition to the DG component near 
the current sheet. This variation of Bjv may result from an 
eigenmode of the current sheet surface waves [Dai, 2009; 
Dai et al., 2011]. The DG Bjv is negative in the northward 
plasma flow and positive in the southward flow, consistent 
with the prediction that the bidirectional plasma flows are 
accelerated by a JxBjv slingshot effect resulting from recon¬ 
nection. The observed bidirectional reconnection outflows 
indicate that the spacecraft moved from the northward side 
to the southward side of an active reconnection X-line. 

Figure 2g exhibits a significant deviation of the mea¬ 
sured electric field Em from —Vi x B from 23:24:44UT to 
23:24:58UT during the current sheet crossing. This devia¬ 
tion illustrates the violation of the ion frozen-in condition 
as a support of the reconnection electric field Em, exper¬ 
imentally defining the so-called ion diffusion region. From 
23:25:20UT to 23:25:44UT is the edge of the diffusion re¬ 
gion and far from the current sheet center. The ion diffu¬ 
sion region in the current sheet is marked by a pink rect¬ 
angle. Within the diffusion region, a remarkable normal 
electric field Ejv as large as 10 mV/m was observed. The 
out-of-plane magnetic field Bm displayed one bump around 
the current sheet center (B^ ~ 0); it was about 10 nT 
larger than the average guide field outside the current sheet. 
The location and waveforms of Ejv and Bm are completely 
consistent with those of the Hall magnetic fields and elec¬ 
tric fields in the asymmetric reconnection reported in early 
THEMIS observations [Mozer et al., 2008; Bonnell et al., 
2008]. Gontrary to the situation of a quadrupole Hall Bm 
and a bipolar Hall Eiv in a symmetric reconnection, there 
is only one bump in the Hall field in the asymmetric recon¬ 
nection. The observed asymmetric Hall electric fields and 
magnetic fields indicate the operation of collisionless recon¬ 
nection. 

The length of the diffusion region along the normal direc¬ 
tion is the crossing time AT multiplied by the magnetopause 
velocity v„ relative to the spacecraft. We estimate v„ by as¬ 
suming that the magnetopause had a constant tangential 
electric field as it moved in the normal direction at constant 
speed. The estimated v„ is [Em(1)-Em(2)]/[Bl[1]-Bl[2]] 
[Mozer et al., 2002]. Here the fields EM(l),EM(2),Bi[l] and 
Bl[2] are measured at the spacecraft frame at times 1 and 
2 upstream and downstream of the current sheet, respec¬ 
tively. As shown in Figure 2d and 2g, Em(1) is ~0.5mV/m 
and Bz,[l] is~ -45 nT from 23:24:20UT to 23:24:40UT be¬ 
fore the crossing. Immediately after the crossing Em (2) 
is ~2.5mV/m and Bi,[2] is~70 nT from 23:25:02UT to 
23:25:14UT. The magnetopause velocity V„ is ~ 17.4 km/s. 
The length of the diffusion region along the normal direction 
is 264.5 km, about 7.3 magnetosheath ion skin depths or 1.2 
magnetosphere ion skin depths (c/ujpi). This scale size is 
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consistent with the values of the diffusion region in previ¬ 
ously reported examples of collisionless reconnection [Mozer 
et al, 2002; Vaivads et al., 2004; Wygant et al., 2005]. The 
estimated tangential reconnection electric field near the re¬ 
versal of Bl is ~0.5-1.3mV in the frame co-moving with the 
magnetopause current sheet. 

The departure of the pressure tensor from cylindrical 
symmetry about the local magnetic field direction can be 
measured by a nongyrotropy index [Aunai et al, 2013].The 
nongyrotropy and similar agyrotropy indexes have been suc¬ 
cessfully applied to characterize non-gyrotropic electrons in 
the electron diffusion region [Scudder et al, 2012; Tang 
et al., 2013; Aunai et al, 2013]. The comparison between 
nongyrotropy and non-ideal electric held, however, has not 
been done yet in these studies. In our event. Figure 2h shows 
an intensified layer of ion nongyrotropy within the diffusion 
region, indicating strongly non-gyrotropical ions when the 
ions are demagnetized. The ion nongyrotropy is as large as 
0.3 in the diffusion region, significantly larger than the aver¬ 
age value (~0.05) outside the diffusion region. Both the off- 
diagonal pressure component and the difference between the 
diagonal component can contribute to the nongyrotropy in¬ 
dex [Aunai et al., 2013]. We hnd that the contribution from 
the off-diagonal pressure component is dominant (>90%) in 
the diffusion region. 

The non-gyrotropic pressure can lead to the violation of 
the ion (or electron) frozen-in condition via the off-diagonal 
terms in the pressure tensor. As similar to the situation of 
non-gyrotropic electron pressure [Vasyliunas, 1975; Hesse 
et al., 1999, 2011], the gradient of the off-diagonal ion 
pressure components can give rise to a non-ideal reconnec¬ 
tion electric field in the diffusion region, (E-|-Vi x B)m ~ 
(dPNM/dN + dPML/dL)/niqi. Figure 2i presents the ion 
pressure tensor components in the boundary normal coor¬ 
dinate system. Spin-resolution ion moments data has been 
extensively used in the Walen test (or tangential compo¬ 
nent test) [Sonnerup et al., 1981]. Ion moments data are 
generally considered OK when the spacecraft takes mul¬ 
tiple samples in the reconnection current layer [Hoshino 
et al., 1998; Nagai et al., 1998; 0ieroset et al., 2001; Mozer 
et al., 2002; Wygant et al, 2005; Zhou et al., 2009; Aunai 
et al., 2011]. From 23:24:46UT to 23:24:54UT, THEMIS 
observed a signihcant non-ideal reconnection electric field 
(E-I-Vi X B)m of 3-4.5 mV/m (Figure 2g) in the center 
of the diffusion region associated with a gradient of the 
ion off-diagonal pressure component (Figure 2i). The den¬ 
sity Hi (Figure 2b) was ~6.2cm“^ at 23:24:50UT, AN from 
23:24:46UT to 23:24:54UT is ~-129km, Pnm (Figure 2i) de¬ 
creases from 2500eVcm“® to near 0, APjvm ~ -2500eVcm“^. 
With all these numbers, our estimate of APnm /^Nuiqi is 
-|-3.2mV/m. APml is ~ -lOOOeVcm”^ in Figure 2i. As¬ 
suming AL~ AN, APML/ALmqi is ~ 1.3mV/m. The as¬ 
sumption of AL~ AN corresponds to an ion diffusion region 
extending 16 ion skin depth (roughly from 23:24:46UT to 
23:25:20UT) northward in the L direction. Such L-extent of 
ion diffusion region is consistent with reconnection models. 
The gradient of the off-diagonal ion pressure terms agrees 
well with the non-ideal reconnection electric held, indicating 
that the non-gyrotropic pressure effect mainly contributed 
to breaking the ion frozen-in condition. Ion inertial (acceler¬ 
ation) effects and, in principle, the anomalous collision effect 
can also contribute to reconnection electric fields. In the dif¬ 
fusion region, Vjv~-70km/s and AVm~130 km/s (Figure 
2a). The ion inertial term scales as {rm/qi)Y nAVm/AT i 
(<0.3mV/m), much less than the non-gyrotropic pressure 
effect. The anomalous collision term is expected to be very 
small, because the observed reconnection should be collision¬ 
less as implied by current sheet thickness (ion skin depth) 


and the Hall helds. The ion pressure tensor can be decom¬ 
posed into a gyrotropic part and a non-gyrotropic part. As 
shown in hgure 2j, the non-gyrotropic part of the pressure 
tensor mostly contributes to the off-diagonal components. 
Based on the analysis of all possible effects in the ion mo¬ 
mentum equation, we demonstrate that the non-gyrotropic 
pressure effect is primarily responsible for breaking the ion 
frozen-in condition in this event. 

The ion momentum equation can be checked term by term 
in the normal direction as well. (E-|-Vi x B)Ar is about 7-10 
mV/m near the reversal of Bi^. The ion inertial term in 
the normal direction is on the order of 0.1 mV/m. The 
gradient of the pressure term in the normal direction is 
{dPNN/dN + dPLNldL)lniqi. The APjvjv ~ -6000eVcm“® 
is the dominant term, contributing to 7.8 mV/m in the 
ion momentum equation. Within a difference of 20%, the 
pressure gradient term is approximately consistent with the 
(E-|-Vi X B)]v. This independent verification of the ion mo¬ 
mentum equation in the normal direction involves the quan¬ 
titative comparison between three independent data sets, 
indicating that THEMIS measurements were reasonably re¬ 
liable in the ion diffusion region. 

As implied from the the two-ffuid momentum equation 
(1), electrons are expected to exhibit similar dynamics as 
ions except that the non-ideal electron terms operate on the 
smaller electron scale in the diffusion region. Theoretical 
studies show that non-gyrotropic pressure effect unfreeze 
electron fluid near the reconnection X-line [Hesse et al., 
1999; Kuznetsova et al., 2001]. Encompassing the obser¬ 
vation of ion non-gyrotropic pressure, an emerging scenario 
is that the non-gyrotropic pressure effect applies universally 
to unfrozen ions at large spatial scale and electrons at small 
spatial scale. Near the X-line, the reconnection electric field 
should be equally expressed in terms of ion or electron non- 
gyrotropic pressure [Hes.se et al, 2011]. 

Figure 3 shows the ion velocity distributions perpendic¬ 
ular to the magnetic field near and in the diffusion region. 
Panel a shows a typical gyrotropic ion distribution in the 
magnetosheath outside the diffusion region. The distribu¬ 
tions in panel b and c are centered around Bi,=-20nT and 
B_L=-|-40nT, where the magnetic held are found to deviate 
from the spin-averaged values by ~30% and ~20% during 
one spin period, respectively. The velocity distributions in 
panels b and c are in the clearly identihed ion diffusion re¬ 
gion and characterized by non-gyrotropic bulges in the core 
component. According to the dehnition of pressure tensor 
P = mf dv^(v — u)(v — u)/(x, v,t), where u is the bulk 
how velocity, such non-gyrotropic bulges are expected to 
contribute to non-zero integral of (vm — UM)(viv — uw) and 
thus large off-diagonal pressure components. Two bulges are 
centered at (50,15) km/s and (-82,-50) km/s in the interval 
c. The velocity spread of these bulges is ~100km/s, about 
half of the thermal speed in the diffusion region. The ob¬ 
served non-gyrotropic distribution is particularly important 
for magnetic held line reconhguration in the diffusion region 
as it can drive a current sheet instability that induces the 
normal magnetic held component ]Mofsc/imann and Glass- 
meier, 1998a]. In addition, the non-gyrotropic distributions 
are expected to drive plasma waves that tend to stabilize the 
plasma [Motschmann and Glassmeier, 1998b]. The strongly 
non-Maxwellian distribution observed in the diffusion region 
has two implications: 1) the effective ion-electron collision 
is insufficient to lead to Maxwellian distributions while ions 
are demagnetized, and 2) the ion demagnetization process 
is much faster than ion thermalization and thus unlikely to 
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be explained by effective ion-electron collisions. The non- 
gyrotropic bulges in panels b and c resemble those reported 
by GEOTAIL spacecraft in reconnection [Hoshino et al, 
1998; Nagai et al., 1998; Wygant et al, 2005]. Numerical 
simulations suggest that non-gyrotropic bulges are caused 
by the mixing of already-accelerated meandering ions with 
ions just convected into the vicinity of the X-line geometry 
[Hoshino et al, 1998; Hoshino, 1998]. Such a mixing process 
is inferred to also occur in the asymmetric reconnection at 
the magnetopause from observations from CLUSTER [Lee 
et al., 2014]. Ions from distinct sources form phase-bunched 
bulges in the distribution. This scenario of forming non- 
gyrotropic bulges is supported by the difference between 
the energy of the non-gyrotropic bulges shown in panels b 
and c. The nongyrotropy in the warmer ion components in 
panel d resemble those in THEMIS observations outside the 
current sheet center but near the reconnection site [Zhou 
et al, 2009]. These non-gyrotropic warmer ions are formed 
by a cucumber-type trajectory that has the neutral sheet¬ 
crossing and the non-crossing segments near an X-line ge¬ 
ometry [Zhou et al, 2009]. 


3. Summary and Conclusions 

Ion thermal pressure effects have been neglected and the 
ion inertial effect is difficult to identity in the generalized 
Ohm’s law. Therefore, the individual ion momentum equa¬ 
tion without approximation is the best choice for investi¬ 
gating the ion-scale physics in reconnection. We take this 
approach and present clear evidence that a non-gyrotropical 
ion pressure situation is responsible for breaking the ion 
frozen-in condition in reconnection. For the first time, we 
show that the non-ideal reconnection electric field associated 
with the violation of ion frozen-in condition is consistent 
with the gradient of the off-diagonal ion pressure compo¬ 
nents. Our approach and results would be valuable for the 
MMS (Magnetospheric Multiscale Mission) mission that is 
designed to resolve all scales relevant to magnetic reconnec¬ 
tion. 

Electrons are expected to exhibit similar dynamics as ions 
except that the non-ideal effect kicks in at smaller electron 
scale. An emerging scenario is that non-gyrotropic pressure 
effect universally unfreezes ions at ion scale and electrons 
at smaller scale in the collisionless reconnection. As evident 
from the two-fluid momentum equation, the reconnection 
electric field is equally explainable in terms of ions or elec¬ 
trons. Near the X-line, the reconnection electric field may be 
equally expressed in terms of ion or electron non-gyrotropic 
pressure effect [Cai and Lee, 1997; Hesse et al, 2011]. 
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Figure 1. Schematics of magnetic reconnection in the 
Earth’s magnetopause, a, The magnetic field geometry 
at the dayside magnetopause reconnection viewed in the 
noon-midnight plane. Magnetic field lines near magne¬ 
topause reconnection can be divided into three classes ac¬ 
cording to their topology : (1) interplanetary field lines 
(red) with no magnetic foot on the earth, (2)‘open’ held 
lines (green) with one magnetic foot connected to the 
earth and (3) ‘closed’ held lines (blue) with both mag¬ 
netic feet connected to the Earth. Two branches of sep- 
aratrix surface SI and S2 intersect along a magnetic ‘X- 
line’ directed out of the plane, b, The expanded view 
of the region surrounding the X-line. A portion of the 
plasma outhow came across the ion diffusion region iden- 
tihed by a signihcant deviation of E from — x B. 
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Figure 2. Observations by THEMIS spacecraft E in the 
ion diffusion region. All THEMIS data in Figure 2 are of 
spin-resolution and at the cadence of the particle instru¬ 
ment. a, Three components of the proton flow velocity 
in the LMN coordinate system, showing a flow rever¬ 
sal in the L (North-south direction), b), The plasmas 
density, c, The electron differential energy flux, d-e, 
three component of magnetic field in the LMN coordi¬ 
nate. f,g, The comparison of v; x B with the electric 
fields in N and M, respectively. Three dimensional elec¬ 
tric fields are obtained from the E ■ B=0 assumption, h 
The ion nongyrotropy index [Aunai et al., 2013], char¬ 
acterizing the degree of nongyrotropy in the distribution 
function. The ion velocity distribution function in inter¬ 
vals marked by A,B,C and D are show in Figure 3. i. 
Components of the ion pressure tensor in the boundary 
normal coordinate, j. The non-gyrotropic part of the ion 
pressure tensor. 
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Figure 3. Observations of ion velocity distributions near 
and within the ion diffusion region, a-d, The cut of ve¬ 
locity distribution functions perpendicular to the mag¬ 
netic field measured from THEMIS ESA in the intervals 
marked as ‘A,B,C,D’ respectively in Figure 2. The hor¬ 
izontal axis is the direction of ion convection flow 
and the vertical axis is the direction of BxVix- The 
ion bulk flow velocity has been subtracted from the dis¬ 
tributions. The white horizontal line in the distribution 
represents the flow convection velocity in the spacecraft 
frame. Interval a is far from the diffusion region, b and c 
are within the diffusion region and d is at the edge of the 
current layer and away from the current sheet center. 
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